Localized collapse and revival of coherence in an ultracold Bose gas 
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We study the collapse and revival of coherence induced by dipolar spin waves in a trapped gas 
of S7 Rb atoms. In particular we observe spatially localized collapse and revival of Ramsey fringe 
contrast and show how the pattern of coherence depends on strength of the spin wave excitation. 
We find that the spatial character of the coherence dynamics is incompatible with a simple model 
based only on position-space overlap of wave functions. This phenomenon requires a full phase-space 
description of the atomic spin using a quantum Boltzmann transport equation, which highlights spin 
wave-induced coherent spin currents and the ensuing dynamics they drive. 

PACS numbers: 03.65.Yz, 51.10.+y, 67.85.-d, 75.30.Ds 



Trapped ultracold gases are promising systems for 
compact, high resolution interferometric sensors. They 
have been proposed for a diverse range of applications, 
including atomic clocks 111, [3] , magnetometers Q , polar- 
izability measurements [j], and inertial force sensors 
Many interferometric sensors rely on differential energy 
shifts for measurements. Furthermore, ultracold trapped 
gases all require inhomogeneous external potentials for 
confinement. While this confinement allows for long in- 
terrogation times and high sensitivities, the inherent in- 
homogeneity can have profound impacts on the coher- 
ence of the sample. Thus it is imperative to understand 
atomic coherence in trapped gases in the presence of in- 
homogeneous potential energies if they are to be used as 
interferometric sensors. 

While Bose-Einstein condensates offer the potential for 
enhanced measurement sensitivity via squeezed states, 
it is often equally advantageous to use ultracold non- 
degenerate gases due to their larger populations, faster 
cycle times, and simpler experimental designs [6(. Even 
above degeneracy, however, gases can exhibit quantum 
scattering effects if the de Broglie wavelength of the 
atoms is long compared to the scattering length a. In this 
case, interferometric measurements using non-degenerate 
quantum gases must contend with the identical spin ro- 
tation effect (ISRE) [7]. Exchange symmetry between 
identical particles in a quantum gas leads to interference 
effects: if the interacting particles have different spins, a 
net spin rotation can occur. This spin rotation gives rise 
to a wide array of macroscopic collective behaviors, in- 
cluding spin-state segregation H @] , spin waves [l(J EH , 
spin locking in coupled fluids [12j . and dramatically en- 
hanced coherence times Q. 

Of particular concern for trapped atom-based measure- 
ments is the collapse and revival of coherence, which was 
first noted in 13Hl5j. A recent study of the collapse 



If the trap frequencies are larger than all other fre- 
quency scales, especially the collision rate and differential 
frequency shifts, as in Qjthe ensemble remains largely 
spin-synchronized. (See 16J for a detailed discussion of 
time scales.) This limit obscures the spatial character 
of the spinor evolution, and for experiments outside this 
limit, a more detailed approach is required. 

In this Letter, we observe localized collapse and re- 
vivals of interference fringes and demonstrate how the 
spatial character of the localization changes with excita- 
tion strength. Furthermore, we show that a simple, intu- 
itive model based on the degree of spatial overlap of the 
spin-state wave functions fails to explain the observed in- 
terference patterns, requiring a full numerical treatment 
using a quantum Boltzmann equation for spin transport. 

We interferometrically study the coherent evolution of 
atomic spins in the presence of spin wave excitations. 
The spin system described herein consists of two mag- 
netically trapped hyperfine states of 87 Rb, |1) — \F — 
l,mp — —1) and |2) — \F — 2 1 mp — 1), which form a 
pseudospin-1/2 doublet. A two-photon microwave tran- 
sition couples |1) and |2) and creates an equal coherent 
superposition. We excite spin waves in this superposition 
with an inhomogeneous differential optical potential Udis 
designed to drive the lowest normal spin wave mode: the 
dipole mode 
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and revival of interference in the time domain showed 
that the dynamics could be tuned by altering the dif- 
ferential potential experienced by the two-state system 



The optical potential consists of a diode laser detuned 
by 0.3 nm from the D% cooling transition, and the slight 
difference in the energies of the hyperfine ground states 
gives slightly different detunings for each state, produc- 
ing a differential AC Stark shift. A frequency and ampli- 
tude modulated acousto-optic modulator sweeps the laser 
across the atom cloud and creates a differential potential 
energy that varies linearly with axial position z (see [ill ] ) . 
The gradient of the spin perturbation, dUdis/dz, is con- 
trolled by adjusting the laser power. The Zeeman shift 
from the trapping potential and the mean-field shift also 
contribute to the differential potential, but these sources 
of inhomogeneity can be eliminated with a mutual com- 



pensation scheme by tuning the magnetic field [8[. 

Typical experimental parameters include a magnetic 
field of Bo = 3.05 G at the center of the trap for this 
mean-field compensation, with a peak density of uq = 
2.5 x 10 19 m~ 3 and temperature of 650 nK, about 60% 
above the critical temperature for condensation. The 
cylindrically symmetric magnetic trap is highly elongated 
in the axial direction with an aspect ration of 37:1, cor- 
responding to trap frequencies of 247 and 6.7 Hz. Thus, 
the density may be radially averaged, and all dynamics 
can be treated one-dimensionally. 

Ramsey interferometry allows us to measure the lo- 
cal transverse components of the atomic spin using two 
7r/2-pulses separated by a variable delay time t. Follow- 
ing the interferometer sequence, absorption imaging mea- 
sures the number of atoms Ni returning to the |1) state. 
The microwave coupling is detuned a small amount from 
the |1) — 1 2) transition (S ~ 30 Hz), and as we vary t we 
observe high contrast interference fringes (Fig. [IJ. Al- 
ternately we can omit the final 7r/2-pulse to project the 
superposition onto the |1) state to measure iVi, or re- 
place the final 7r/2-pulse with a it- pulse to measure N%. 
For all measurements, we measure the atomic distribu- 
tions with absorption imaging and bin the images into 
equally spaced axial bins. 

The optical potential excites a dipole spin wave mode 
that causes a linear phase gradient to build up across 
the cloud, as seen by the tilting of the phase fronts in 
Fig.QJa). ISRE-induced spin currents build up and lead 
to longitudinal rotation of the spin. This spin rotation 
manifests as axial separation of the spin states in the 
magnetic trap [Fig. Hfd)], and as the spin wave peaks, 
the phase fronts become disrupted and interference fringe 
visibility collapses. The spins soon rephase, and the spin 
oscillation continues, with a second collapse occurring at 
the next peak in the spin wave cycle. 

There is an overall damping of both the spin wave ex- 
citation and the superposition coherence through elas- 
tic collisions that eventually randomize the atomic spins 
and destroy the coherent superposition. This damping is 
more pronounced in potentials with greater inhomogene- 
ity. There is also longitudinal relaxation due to dipolar 
collisions leading to loss from the |2) state. Lastly the 
spin waves themselves damp by an additional mechanism: 
Landau damping due to coupling to higher order, highly 
damped spin wave modes [10|, LL3j, |l4 1 . 

The spatial pattern of the interference contrast is par- 
ticularly interesting, as an ensemble-wide collapse does 
not occur. There is nearly complete collapse in the cen- 
ter of the atom cloud [Fig. H^b)], but in the wings of the 
distribution, fringe visibility is only moderately affected 
at the peak of the spin wave perturbation [Fig. HJc)]. 
Furthermore, this spatial pattern of interference is not 
identical for all spin wave-induced collapses. Figure[2]dis- 
plays Ramsey fringes for a more weakly driven spin wave 
excitation. When the fringe collapse occurs, it happens 
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FIG. 1. (color online) Ramsey fringes in the presence of 
a dipolar spin wave driven by an optical potential with 
dUdif{/dz ~ ft x 10 Hz/mm. (a) Ni measured after a Ramsey 
sequence, radially averaged and normalized to its maximum 
value, and plotted in false color as a function of t and z. Data 
is interpolated between bins and time steps, (b)-(c) Cross 
sections of Ramsey fringes for the center bin and a bin near 
the edge of the cloud respectively, as shown by the arrows in 
(a), (d) Normalized projections of the of the states of the su- 
perposition N\ (solid line) and N2 (dashed line) at the time 
of maximal component separation and collapse, denoted by 



an arrow m 



(a). 



to a greater degree near the edges of the cloud and not as 
much in the center [Figs.EJc) and (b) respectively]. Note 
that the time scales are different, as spin wave frequency 
has a nonlinear dependence on excitation amplitude in 
this system [Til ] . 

To highlight the need for a complete phase-space anal- 
ysis of the spin dynamics, we construct a simple model 
for fringe visibility F Q based on the spatial overlap of the 
two spin states in the superposition. The local relative 
populations of the superposition, Ni(z) and N 2 {z), are 
projected as described previously [e.g. Fig. QJd)]. If the 
interference contrast were determined only by the local 
density overlap, then it would be given by 



F Q {t,z) 



2 y /N 1 {t,z)N 2 {t,z) _ f/T 
N 1 {t,z) + N 2 (t,z) € ' 



(1) 



where r is an empirically determined time constant that 
characterizes the cloud-wide decoherence due to elastic 
scattering. We set t ~ 350 ms using an average of the de- 
coherence times over the entire distribution. Figure [3ja) 
shows how F a evolves in the presence of the same dipolar 
spin wave perturbation shown in Fig. [TJ while Fig. E|b) 
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FIG. 2. (color online) Ramsey fringes in the presence of 
a more weakly driven dipolar spin wave (dU^jg/dz — h x 
7 Hz/mm), (a) iVi measured after a Ramsey sequence, nor- 
malized to its maximum value, (b)-(c) Cross sections of Ram- 
sey fringes for the center bin and a bin near the edge of the 
cloud respectively, as shown by the arrows in (a), (d) Nor- 
malized projections of Ni (solid line) and N2 (dashed line) 
at the time of maximal component separation and collapse, 
denoted by an arrow in (a). 



shows the actual measured fringe contrast. The mea- 
sured contrast is more sensitive to shot-to-shot magnetic 
field-induced phase noise, and the reduced contrast at 
~ 150 ms on the positive side of the trap (where the \2) 
state density is highest) is caused by atom loss due to 
dipolar relaxation. 

The large component separation means the longitu- 
dinal rotation of the spin is especially pronounced on 
the edges of the cloud during the spin wave peak (i ~ 
180 ms), and thus the density overlap method predicts 
the smallest fringe contrast to be there. In the center 
of the atomic distribution the density overlap is maxi- 
mal, and one would expect higher contrast than on the 
edges [Fig.|3fa)]. This prediction clearly does not explain 
the actual collapse and revival dynamics observed in the 
strong dipole potential, which gives minimal contrast in 
the cloud center during the spin wave peaks [Fig. |3]Jb)]. 
While the local density overlap does give a maximum pos- 
sible value for interference contrast at a given position, 
it provides an incomplete picture. 

We turn to a full phase-space treatment of the spin 
vector to understand the spin dynamics completely. The 
spin evolution can be described with a quantum Boltz- 
mann transport equation for the spin momentum and 
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FIG. 3. (color online) (a) Predicted interference contrast F Q 
determined from population overlap (see text), (b) Measured 
fringe contrast, showing minima in the cloud center at t ~ 180 
and 575 ms, calculated by performing Fourier transforms on 
the Ramsey fringes in Fig. OJa) using a two-fringe running 
window to extract the fringe amplitude. The resulting con- 
trast is filtered to remove residual aliasing at the Ramsey 
frequency. 
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where m is the atomic mass, £/ cx t is the axial harmonic 
trapping potential, and the right hand side represents a 
collisional relaxation term related to the elastic scatter- 
ing rate. The differential potential Udis and mean-field 
coupling g that drive spin waves are contained in the 
cross term x a as follows, 



£1 = (gS + Udisw)/h. 



(3) 



The mean-field coupling constant is g = Airh 2 a/ra, and 
w is a longitudinal unit vector on the Bloch sphere. The 
spatial spin distribution is obtained by integrating a over 
the momentum distribution, S(z) — J adp. 

We numerically solve Eq. [2] using an alternative di- 
rection implicit finite difference method. The transverse 
component a± describes the coherence. Figure 2] shows 
the magnitude of a± and its orientation <j) for a selection 
of evolution times, as well as the Ramsey fringe contrast 
obtained by integrating a± over p to find \S±\. 

The full phase-space evolution of a± sheds light on the 
spatial nature of the coherence as follows. The effect of 
Udis manifests in Fig. [4] as an axial phase gradient (ver- 
tical phase stripes) for t < 100 ms. As the atoms move 
in the confining potential, these phase gradients correlate 
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FIG. 4. (color online) False color images from numerical simulations of the evolution of <tj_ (top) and cf> (middle) in phase space 
under the influence of a dipolar spin wave, using similar parameters as for the data shown in Fig. [T] This initial distribution at 
t — is Gaussian in p and z. At each time step a± is normalized and does not reflect the global relaxation of the spin vector. 
Bottom row: interference fringe contrast (transverse spin distribution S±) obtained by integrating a± over the momentum 
coordinate. The vertical scale for t > 70 ms has been magnified by a factor of 3 to compensate for transverse spin relaxation. 



with momentum and position, and coherent spin currents 
develop on either side of the trap (t ~ 100 — 120 ms). 
The ISRE keeps the phase within each lobe of the spin 
distribution more uniform than the differential potential 
would otherwise suggest. The two lobes have opposite 
phases from each other, however, and where they meet, 
nodes form in the spin distribution. When the counter- 
propagating spin currents meet at the peak of the spin 
wave cycle, exchange scattering rephases the two halves 
of the distribution {t ~ 130 — 140 ms, corresponding to 
maximum spin-state segregation). If the phases of the 
two lobes are significantly different, this rephasing leads 
to a zero in the center of the transverse spin distribution 
and manifests as a localized collapse of coherence. Co- 
herence revives somewhat as the spin wave completes a 
cycle (t > 150 ms) due to the ISRE rephasing the entire 
spin distribution, although elastic scattering destroys the 
superposition on a similar time scale. 

For smaller amplitude spin wave excitations, the nodes 
develop more slowly and do not reach the center by the 
time the spin wave peaks. Instead they are ejected from 
the phase-space profile, and thus the regions of minimal 
contrast do not reach the center of the distribution, as 
in Fig. [5] Conversely, if the spin wave amplitude is in- 
creased further, phase accrues more rapidly, and more 
nodes begin to appear. The presence of contrast minima 
in the wings of the distribution at t ~ 180 ms in Fig.[3Jb) 
may indicate that we are beginning to enter this regime. 
Decoherence from elastic scattering increases as well for 
more strongly driven spin waves, and thus detection of 
more complex spatial features is difficult. 

We have demonstrated that the coherence of trapped 
gas interferometers near quantum degeneracy can exhibit 
complex spatial patterns of collapse and revival. This be- 
havior can be understood by considering the full phase 
space evolution of the spin distribution given by a quan- 



tum Boltzmann equation. Although we cannot tomo- 
graphically reconstruct the full a distribution, the effects 
of quantum dynamics in the momentum coordinate are 
still manifested in the spatiotemporal coherence patterns. 
Clearly care must be taken to avoid significantly inhomo- 
geneous differential potentials when attempting to make 
interferometric measurements. Such inhomogeneous po- 
tentials enhance dephasing, and even worse, can lead to 
complete collapse of coherence. If inhomogeneous po- 
tentials are unavoidable, their effect may be lessened by 
working with tight confining potentials and lower densi- 
ties, as in Q, where the entire atomic distribution tends 
to stay phase synchronized to a greater extent. 
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